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ABSTRACT 

The standard cosmological model based on cold dark matter (CDM) predicts a large number of 
subhalos for each galaxy-size halo. It is well known that matching the subhalos to the observed 
properties of luminous satellites of galaxies in the local universe poses a significant challenge to our 
understanding of the astrophysics of galaxy formation. We show that the cosmic evolution and 
host mass dependence of the luminosity function of satellites provides a powerful new diagnostic to 
disentangle astrophysical effects from variations in the underlying dark matter mass function. We 
illustrate this by comparing the results of recent observations of satellites out to z = 0.8 based on 
Hubble Space Telescope images with the predictions of three different sets of state-of-the art semi- 
analytic models with underlying CDM power spectra and one semi-analytic model with an underlying 
Warm Dark Matter (WDM) power spectrum. We find that even though CDM models provide a 
reasonable fit to the local luminosity function of satellites around galaxies comparable or slightly 
larger than the Milky Way, they do not reproduce the data as well for different redshift and host 
galaxy stellar mass. This tension indicates that further improvements are likely to be needed in the 
description of star formation if the models are to be reconciled with the data. The WDM model 
matches the observed mass dependence and redshift evolution of satellite galaxies more closely than 
any of the CDM models, indicating that a modification of the underlying power spectrum may offer 
an alternative solution to this tension. We conclude by presenting predictions for the color magnitude 
relation of satellite galaxies to demonstrate how future observations will be able to further distinguish 
between these models and help constrain baryonic and non-baryonic physics. 
Subject headings: galaxies: halos — galaxies: structure — galaxies: dwarf — dark matter 



1. INTRODUCTION 

The standard model based on dark energy and Cold 
Dark Matter (CDM) is successful at describing the 
large-scale ( > Mpc ) structure of the universe (e.g. 
iKomatsu et al.|[2011[ ). In this standard model, galaxies 
and clusters of gal axies grow hierarchi cally within dark 
matter halos (e.g., [Spring el et al]l200fih . 

At the scales typical of galaxies (~kpc), the model 
faces a number of observational challenges. A major 
source of tension between observation and the stan- 
dard model is the s o-called "missing s atellite problem" 
(jKlvpin et all 119991: iMoore et allll999t) . At large scales, 
a simple ranked matching of simulated subhalos to ob- 
served luminous structure predicts clustering behavior 
consi stent with that of the observed luminous s truc- 
ture (|Conrov fc Wechslerll2009HBehroozi et~aT]|2010D . At 
smaller scales however, this simple matching scheme 
breaks down. For instance, CDM simulations predict 
that the Milky Way (MW) should have approximately 
10 subhalos with rotation velocities greater than Fornax, 
while only 3 satellites that are more luminous are ob- 
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serve d (jBovlan-Kolchin et al.l 120121 : IStrigari fc Wechsleri 
.2012-!). The discrepancy between the predicted number 
of subhalos and the observed number of luminous satel- 
lites becomes more d ramatic at even lower halo masses 
(Stri gari et ail l2007t ) . While the measurement of low 
luminosity satellites may be biased in the M ilky Way 
due to disk obscuration ([Tollerud et al.l [20081. the lack 
of low mass galax i es ext ends to the field as well, where 
iPapastergis et al.l (|2011| ) measured a factor of 8 fewer 
galaxies with velocity widths of ^50 km/s than predicted 
by current ACDM simulations. Furthermore, recent sur- 
veys have found that there is not a si gnificant populatio n 
of optically dark, gas rich galaxies, (iDovle et al.l r2005). 
indicating that if these very low mass halos exist, they 
do not contain significant amounts of gas or stars. 

The discrepancy may be due to a variety of factors. 
The complex physics of star formation in low mass ha- 
los makes matching simulated halos to observed satellite 
galaxies difficult. Significant progress has been made in 
understanding the complex p rocesses which affect star 
formation in low mass halos (|Kravtsovl 120101) . In par- 
ticular, numerous studies have focused on studying how 
both sub-galactic effects such as supernovae feedback and 
stellar winds, and super-galactic effects such as UV heat- 
ing from reionization, and tidal and ram pressure strip- 
ping by the central galaxy can suppress star formation 
and thereby pro duce the present day satellite luminos- 
ity function (e.g I Bullock et al.ll2q00t [Benson et al.ll2002l; 
Somervillel 120021: iKravtsov et al.ll2004 iKaufmann et al.l 



2008|; LMaccio 20101: lSpringelll201ol)" 



In addition to suppressing star formation, thereby 
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making subhalos undetectable, several studies have ex- 
amined whether supernovae feedback can significantly 
flatten the central regions of dark matter density pro- 
files. This would potentially explain the discrep- 
ancy between the predicted and observed central ve- 
locities of ACDM subhalo s around the Milky Way 
(jBovlan-Kolchin et al.l [2012), by lowering the subhalo 
central velocities, thereby making them appear less 
massive. Simulations with baryons have produced 
conflicting predictions a s to whether supern ovae can 
(iGovernato et al.l 120121: iTeyssier et al.l I2013D or can- 
not ( Garrison-Kimmel et al.ll2013l ) significantly alter the 



kinematics in the central regions of dark matter halos. 

As all of these baryonic processes occur over cosmo- 
logical time scales, ideally the predicted effects would be 
compared with observations of the satellite population 
over as much of the history of the universe as possible. 

If the solution to this problem is astrophysical in na- 
ture, i.e. low mass halos exist but do not form stars 
and are therefore undetected by traditional astronom- 
ical observations, one needs to turn to other methods 
to verify their existence. This line of reasoning has 
motivated searches for satellites halos based on proper- 
ties that are independent of their ste llar content, such 
as their gravitational lensing effect (|Mao fc Schneiderl 
19981: 1' Mctcalf & Madau HgOp: iDalal fc Kc-chanekl I2002T 
Amara et al.l 120061: iKeeton fc Moustakasl 120091: iTreul 
20iri iVeeetti et al.l 120121 e.g..). or on their expected 
DM annihilation signal dKuhlen et al.ll2008t iPorter et al.l 
1201 It iStrigaril l20ia lAbazajian fc Kaplinghatl I2012T ). or 
their influence of t idal streams in the Milky Way 
(|Carlberg et al.ll2012T ). Direct detection of dark subhalos 
would be a stunning confirmation of the standard model. 

However, this is not the only possible solution to 
the missing satell ite problem. Wa rm Dark Mat- 
ter (WDM) (e.g. IColombi et al.l 119961 and references 
therein) is an interesting alternative to CDM, poten- 
tially offering an elegant astroparticle solu tion to the 
missing satellite problem ()Lovell et al.ll2012l ). In WDM 
scenarios, small-scale structure is suppressed relative to 
CDM. There are a variety of mechanisms which can 
achieve this, for instance either by reducing the mass 
of the dark matter particle for a thermal relic (e.g 
iSteffenl l2tM ) or by introducing non-thermal particles 
such as sterile neutr inos produced from oscillations of 
activ e neutrinos (e.g.lOlive fc Turned ll982l:IShi fc Fulled 
119991 lAbazaiian et all 120011 : iDolgov fc Hansen! I2002T ). 
The suppression of small-scale structure in turn af- 
fects the full merger history of a halo and its sub- 
halos, which recent cos mological simulati o ns have be- 
gun t o study in detail (iMenci et al.l 120121: lLovell et al.1 
2012t iKamada et al.l 120131; iPolisenskv fc Ricottil 120111: 



Bode et al.ll2001fc iGottloeber et al.ll2010D . 

Clearly the astrophysical and astroparticle solutions to 
the missing satellite problem are not mutually exclusive, 
and progress on both fronts may be needed to reconcile 
the model with the data. In fact, a proper observational 
test of the WDM model predictions for the number of 
satellite galaxies requires an accurate description of the 
astrophysical effects related to star formation. 

In this paper we present a new test of three semi- 
analytic galaxy formation models, with one model im- 
plemented for both a Warm and Cold Dark Matter sub- 
structure mass function, in order to show how variations 



in the halo mass function produce different predictions 
from variations in baryonic physics. We compare the pre- 
dictions from these models with the observed abundance 
of satellite galaxies as a function of host galaxy mass and 
cosmic time. 

This test is made possible by two recent developments. 
On the observational side, the implementation of power- 
ful statistical tools to detect and count satellites in deep 
archival Hubble Space Telescope images provides data to 
comp are with model predictions (|Nierenberg et al1l2011l 
|2012| ). On the theoretical side, advances in computa- 
tional methods have allowed for cosmological simulations 
with unprecedented volume and resolution. In this work 
we focus on predictions from four independent cosmologi- 
cal simulations, which have semi-analytic models applied 
to dark matter merging trees, which we describe in more 
detail below. 

The paper is organized as follows: In Sj^l we describe 
our observations of the satellite luminosity function. In 
$3] we summarize the key aspects of the theoretical mod- 
els. In <21 we compare the observations with the theo- 
retical predictions. In fj5] we compare a few of the main 
properties of the models and present new predictions for 
the distribution of colors of satellite galaxies across cos- 
mic time. In Sj6]we compare the model results with other 
observations taken from the literature. Finally in <J7]we 
conclude with a discussion and summary of the results. 

2. OBSERVATIONS 

Measurements of t h e sate llite population were obtained 
bv iNierenberg et al.l (|2012f ). using ACS 1814 imaging of 
the COSMOS field, where a full description of the anal- 
ysis is given. 

For convenience of the reader, the detection technique 
is illustrated in Figure [2] First, The depth and high 
angular resolution of Hubble images was used to iden- 
tify satellites as much as a thousand times fainter than 
their host galaxies, and out to redshift z = 0.8. After 
adding simulated satellites near bright central galaxies, 
and testing object recovery and accuracy, the analysis 
was restricted to objects with I814< 25 mag AB, in or- 
der to ensure completeness. 

Usin g a single band of photometry, INierenberg et al.1 
(|2012h model the number density of objects near the 
chosen host galaxies as a combination of a uniform den- 
sity of background/foreground objects, in addition to 
a population of satellite galaxies which have a number 
density which increase s radially a s a power law near 
the host galaxies (e g I Chenl 120081: IWatson et all 120121 : 
INierenberg et al.ll2011l l2012f) . The power of the radial 
distribution as well as the number of satellites within a 
fixed magnitude offset from the host magnitude are si- 
multaneously inferred to yield a cumulative luminosity 
f unction. 

INierenberg et al.1 (|2012h showed that the inferred satel- 
lite luminosity function for low mass host galaxies is fully 
consistent with the luminosity function of Milky Way 
satellites, and the s atellites of other similar mass galax- 
ies at low redshift (IStrigari fc Wechsledl2012l : iLiu et all 
120111; iGuo et al.l l2011al).~ Taking into account the differ- 
ences between global background subtraction rather than 
the locally estimated background used in our measure- 
ment, which can cause a factor of two increase in the 
inferred number of satellites due to correlated structure 
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around m assive host galaxies, ou r results are also consis- 
tent with lWang fc W hite (20121) at low redshifts and for 
the most massive satellites (\og\L „/L h ] > —1), our re- 
sults are consistent with those of ([Newman et al.|[2012f) 
at higher redshifts. 

As we will show below, these new constraints on host 
mass and redshift dependence of the satellite luminosity 
function provide distinguishing power between the phys- 
ical effects of varying the halo mass function and star 
formation physics on the satellite luminosity function. 

3. THEORETICAL MODELS 

We study three different semi-analytic models imple- 
mented in three different cosmological Cold Dark Mat- 
ter simulations. In addition, we apply one of the semi- 
analytic models to a Warm Dark Matter cosmology. This 
section provides a description of both the Cold and Warm 
Dark Matter simulations and their models for galaxy for- 
mation. 

3.1. Cold Dark Matter Models 

When studying the satellite luminosity function, it is 
extremely useful to compare multiple cosmological simu- 
lations along with different star formation parameters in 
order to understand the range of theoretical predictions 
given different assumptions. In this work we focus on 
three Cold Dark Matter cosmological simulations, each 
with its own semi-analytic model for galaxy. The simu- 
lation parameters governing star formation and the un- 
derlying dark matter mass function were fully specified 
prior to performing this test and are the same as in pre- 
vious pape rs, in which a more complete description can 
be found (jMenci et all l2012t iLu et al.1 l2012t iGuo et al.1 
I2011bl) . Below we briefly summarize the relevant aspects 
of these simulations. 

3.1.1. Menci et at 2012 

The semi-analyt ic models are described in detail by 
I Menci et al.l ((2012). In brief, the backbone of the models 
are computations of dark matter merging histories that 
can resolve halos down to masses of 10 5 M©, allowing for 
an extremely accurate characterization of the survival 
and merger histories of satellites. This is essential for 
studying the faint end of the satellite luminosity func- 
tion. The models predict mass distributions of halos, 
and the merging rates of DM halos and subhalos based 
on an input power spectrum, and generate luminosities 
for these halos based on a set of star formation prescrip- 
tions. 

These galaxy formation prescriptions connect the com- 
plex baryon physics of star formation to the dynamical 
evolution of the DM haloes. For each subhalo hosting a 
galaxy, the model predicts radiative gas cooling, ensuing 
star formation and associated feedback from Supernovae 
(SNae) events. Cooled gas settles into a rotationally sup- 
ported disk with radius and rotational velocity related 
to the DM mass of the subhalo. This gas gradually con- 
denses into stars at a r ate consistent wi th the observed 
Kennicut-Schmidt law (K ennicuttl fl998) . SNae return 
part of the cooled gas to the hot gas phase at the virial 
temperature of the halo. Star formation is also triggered 
by galaxy-galaxy interaction-driven starbursts, which in- 
duce gas accretion onto the central supermassive Black 
Holes. 



It is important to note that the set of parameters used 
in these models is much smaller than the number of ob- 
servables they are simultaneously consistent with. In 
brief, the star formation parameters are tuned to simul- 
taneously match the field luminosity function in multi- 
le bands from redshifts to 6, the Tully-Fisher relation 
Tullv fc Fish"enll977D . the stellar mass function of field 
galaxies from between redshifts and 4, the stellar mass- 
star formation relation to a redshift of 2, the colors and 
color magnitude relations of both field and cluster galax- 
ie s in addition to the AGN luminosity function. 

iNierenberg et al.1 (|2012[ ) required host galaxies to be 
outside of R200 of any galaxy with a higher stellar 
mass, where R200 w as estimated using the formula from 
iDutton et al.l (|2010t ). Ideally this selection would be du- 
plicated when comparing with the simulatio n, however 
this is not possible for the lMenci et al.1 (|2012l ) models, as 
they do not contain spatial information. This should not 
significantly affect our results, as comparisons with the 
Millennium and Bol shoi N-Body simu l ation s (described 
below) show that the INierenberg et al.l (|2012l ) central iso- 
lation criteria are efficient at returning central host galax- 
ies (~ 90%), and that the use of the matching selection 
criteria in these cases did not significantly alter the in- 
ferred predicted number of satellites per host, relative to 
simply using the simulation selection of central galaxies. 

3.1.2. Guo et al. 2011 



The semi-analytic model of lGuo et al.l (|2011bh is based 
on two very large dark matter simulations, the Mil- 
lennium Simulation (MS; ISpringel et al.l l2005fl and the 
Mille nnium-II Simulation (MS-II; iBovlan-Kolchin et al.l 
I2009D . The box size of the MS is SOO/i^Mpc and its 
merger trees are complete for subhaloes above a mass 
limit of 1.7 x lQ w h- l M Q . The MS-II follows a cube 
of side 100/i _1 Mpc, but with 125 times better mass 
resolution than the MS (subhalo masses greater than 
1.4 x 10 8 /i- 1 M Q ). Both si mulations adopt the WMAP1- 
based ACDM cosmology (jSpergel et al.ll2003l) with pa- 
rameters h — 0.73, Q m = 0.25,17a = 0.75, n = 1 and 
as = 0-9. D u e to the lower resolution of the MS, 
IWang fc White! (|2012l ) found that the luminosity func- 
tions of satellites flattens for satellites with M r > —18. 
For brighter satellites with M r < —18, the simulations 
are consistent. Since satellite luminosity functions will 
be measured to about three-orders of magnitude fainter 
than central primaries in this paper, we will focus on the 
semi-analytic model implemented on MS-II h ereafter. 

In ge neral, the galaxy evolution mode l of I Guo et al. | 
201 lbf) is based on those developed in ISpringel et aTl 
2005h : ICroton et al.l (|2006| ): lDe Lucia fc Blaizotl (|2007| ). 
The model includes a few main modifications such as 
the different definition of satellite galaxies, the gradual 
stripping and disruption of satellites, a mass-dependent 
model of supernova feedback, a modified model for reion- 
ization and a more realistic treatment of the growth of 
stellar and gaseous disks . Free parameters of these mod- 
els in (jGuo et al.ll201~lll were determined to give close 
predictions to the abundance and clustering of low red- 
shift galaxies, as inferred from SDSS, and are functions 
of their stellar mass, luminosity and color. 

There are two types of satellites in the simulation: i) 
those with an associated dark matter subhalo (type-1) 
and ii) those whose dark subhalo has fallen below the 



4 



Nierenberg et al. 




Ax (R host ) 



Fig. 1. — Demonstration of our image processing technique. The upper panels show the original Hubble Space Telescope images centered 
on main galaxies, while lower panels show the host-subtracted images. Red circles show objects we detect near the hosts. While we study 
the satellite population out to much larger radii (about 5 times larger than the images), the host subtraction in the central region is essential 
for allowing an accurate characterization of the satellite spatial distribution, which in turn allows us to isolate the properties of the satellite 
population. Some objects that are too close to the edges, or which are below our detection threshold are not circled. Objects very close to 
the centers of the host galaxies are excluded as the host subtraction becomes inaccurate in the inner regions. 



resolution limit of the simulation (type-2). For the latter, 
the position and velocity of the orphan galaxy is given by 
those of their most bound particle. Type-2 satellites are 
removed from the galaxy catalogues when one of these 
two conditions is fulfilled: 1) the time passed from the 
disruption of the subhalo is longer than their estimated 
dynamical-friction timescale, or 2) the integrated tidal 
forces from the host halo exceed the binding energy of 
the galaxy. 

Here we use the data d ownloaded from 
Ihttp: / /www.mpagarching.mpg . de / mi llennium! for 
registered users. We project the simulation box in three 
orthogonal directions (along their x, y and z axes). In 
each projection we assign each galaxy a redshift based 
on its "line-of-sight" distance and peculiar velocity We 
select isolated primaries using criteria which are directl y 
analogous to those used in iNierenberg et al.1 ([2012). 
Satellites arc defined to be all companion galaxies 
whose distances to the isolated primaries are smaller 
than the halo virial radius (R,2oo)- To directly compare 
results based on COSMOS, the real R v i r provided in the 



database is ignored, and i nstead we used the empirical 
formula of iDutton et al.l ((2010) to estimate M v j r and 
R200 from stellar masses of galaxies. 

3.1.3. Lu et al. 2012 

The baryonic processes impl emented i n thi s sem i- 
analytic model are described in iLu et al.l (|2011l l2012t) . 
Different from the previously published versions, this ver- 
sion of the model is applied on a set of halo merger trees 
extracted from a large cosmological iV-bod y simulation, 
the Bolshoi simulation ([Klvpin et al.|[20Tll) . which has a 
box size 250/i -1 Mpc on a side. The s imulation adopted a 
cosmology favored b y WMAP7 data flJarosik et al.ll2011f) 
and W MAP5 data (jPunklev et al.ll2009t iKomatsu et al.1 
l2009h with parameters with fi mi0 = 0.27, Oa,o = 0.73, 
n h o = 0.044, h = 0.70, n = 0.95 and ct 8 = 6.82. The 
mass resolution of the simulation is 1.35 x lO 8 ft. _1 M0, 
which allows us to track halos and subhalos with mass 
~ 7 x lO 9 ^" 1 !^©. Dark matter halos a nd subhalos are 
ident ified with the Rockstar halo finder (jBehroozi et al.1 
120121 ) based on adaptive hierarchical refinement in phase- 



Evolution of satellites in cold and warm dark matter 



5 



space. 

As with the two other models in this work, the semi- 
analytic model follows the dark matter merger tree and 
calculates the rates of gas cooling, star formation, outflow 
induced by star formation feedback, and galaxy-galaxy 
mergers. The kinematics of satellite galaxies is followed 
by using subhalo information from the simulation when- 
ever the subhalo is resolved. When a halo becomes a 
subhalo, we instantaneously strip the hot gas associated 
with the halo, while the stellar mass and cold gas mass 
remain intact. 

When the subhalo is no longer resolved in the 
simulation, the model applies dynamical friction 
(jBinnev fc Tremainelll987l ) to estimate when the satellite 
has merged into the central galaxy, and it assumes that 
the tidal stripping is strong enough to also strip the stel- 
lar mass and cold gas mass. At this point, the entire cold 
gas disk is stripped and is mixed into the hot gas of the 
host primary halo. Starting from that time, a fraction 
of stellar mass is tidally stripped per orbital timescale. 
The efficiency of the tidal stripping is controlled by a pa- 
rameter, which is tuned to yield 30% of the stellar mass 
is stripped in every orbital timescale to match the con- 
ditional stellar mass function of local galaxies. 

Other parameters governing star formation and feed- 
back are tuned using an MCMC optimization to match 
the local galaxy stellar mass function (Moustakas et al 
2013). The model is guaranteed to produce a global 
galaxy stellar mass function which provides the best pos- 
sible match to the data between 10 9 and 10 12 M Q at 
rcdshift zero, within the observational uncertainty and 
given the chosen model parametrization. 

To compare with observation, hosts were selected using 
the sa me isolation criteria as used in iNierenberg et aLl 
(2012), with hosts required to be not within R200 of 
a host wit h higher stellar mas s, where R200 is esti- 
mated from iDutton et al.1 (|2010f ) , using the relationship 
for early-type hosts. The satellite luminosity function 
was measured within this region, rather within R v ; r as 
given by the simulation. 

3.2. Menci Warm Dark Matter Model 

Warm Dark Matter has been proposed as a means of 
suppressing the satellite luminosity function by reducing 
the number of low mass dark matter halos, suppressing 
power with respect to CDM below a certain cutoff-scale . 
In this work, we use WDM merger trees from lMenci et al.l 
(HoH) which are based on a cutoff scale of ~ 1 Mpc, 
corresponding to a ~0.75 keV/c 2 thermal relic, which 
suppresses the power spectrum on sub-Mpc scales and 
has behavior equivalent to that of a ~ 3.4 keV/c 2 sterile 
neutrino. We chose this mass to be low enough to affect 
satellite-scale structure while still agreeing with limits 
from observations of large-scale and local group structure 
which constrain the particle mass to be larger than 0.6 
keV/c 2 for a thermal r elic and 2.5 keV/c 2 for a sterile 
neutrino DM particle (IViel et al.l l2009t iBovarskv et al.l 
12001 iPolisenskv fc Ricottill2011lL ^ 

The fiducial mass function from I Menci et al.l ()2012l ) 
was b ased on an Exten ded Press Schechter formalism 
(EPS. iBonTeTaLl [1991 . In our predictions, we con- 
sider the effects of a complete suppression of progeni- 
tors with masses below that corresponding to the free- 
streaming scale; this maximizes the possible effects of 



different window functions and collapse thresho lds in the 
WDM merging trees (see iBenson et al.l 120131 ) . We in- 
clude the possible range of model predictions under dif- 
ferent assumptions for the building up of merging trees 
in the uncertainty regions i n Figure [2) 

The mass function from I Menci et al.l ((2012) is based 
on on the Exten ded Press Schechter formalism (EPS, 
iBond et aLlll991|) , which is modified to take into account 
the suppression of structure below the free-streamin g 
scale (see, e.g. iMenci et al.l 120121 : IBenson et al.l l2013h . 
The uncertainties associated with this modification are 
reflected in the width of the prediction for the final WDM 
satellite luminosity function. 

Our goal in this comparison was to explore how the ef- 
fects of varying the power spectrum compare with the ef- 
fects of varying semi-analytic star formation parameters. 
To achieve this, the same semi-analytic model of star 
formation is used for the WDM models as was applied 
to the CDM merging trees in the Menci CDM model. 
The CDM semi-analytic model was selected to provide a 
good match to the field color distribution. IMenci et al.l 
(|2012[ ) show that this model provides a good fit to the 
field luminosity function when applied to a WDM power 
spectrum. 

4. RESULTS 

The improvements in observations and the implemen- 
tation of semi-analytic models in cosmological scale sim- 
ulations allows a comparison between the observed and 
predicted number of satellites as a function of host galaxy 
mass and cosmic time for the first time, thus allowing for 
significantly more discriminatory power than tests based 
only on the MW or the local volume. 

The results of the comparison are presented in Fig- 
ure [2j where we plot the observed number of satellites 
as a function of the ratio between host and satellite lu- 
minosities in SDSS-r, along with the WDM and CDM 
model predictions. Among the CDM model predictions, 
two main trends are evident. First, no one CDM predic- 
tion precisely matches the observation at all redshift and 
stellar mass intervals. This highlights the importance of 
using data from a range of redshift and stellar mass when 
tuning the parameters of semi-analytic models. All three 
models show qualitatively similar behavior with a strong 
dependence on the number of satellites per host on the 
host galaxy mass. Taking into account the covariance be- 
tween data points (calculated by bootstrap resampling), 
the generalized chi-squared between the models and the 
data is 529, 1960 and 101 for the Menci, Guo and Lu 
models respectively for 41 degrees of freedom. 

The Menci WDM model shows distinct behavior in 
comparison with the CDM models. Notably, it predicts 
weaker host mass dependence and less redshift evolution 
than predicted by any of the CDM models. Of the four 
models it provides the best agreement with the data with 
a chi-squared of 59. 

One of the most interesting results of the above com- 
parison is the significant difference in the predicted satel- 
lite luminosity functions, even among the three CDM 
models. In the following section, we perform a detailed 
comparison of some of the properties of the four models 
in order to explore the cause of these differences. 

5. COMPARISON OF SIMULATION PROPERTIES 
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Fig. 2. — Comparison of CDM and WDM model predictions to observed number of satellites over a range of redshift and host stellar 
masses. The purple, blue and red curves represent predictions from three CDM based simulations with separate semi-analytic models for 
star formation. The gray curve shows the prediction for one WDM simulation with star formation parameters identical to those used in 
Menci, CDM prediction. For the Lu and Guo models, the line widths represent the scatter in the mean value across the simulation boxes. 
For the Menci CDM model, the shaded region represents the range of model predictions possible given the observational uncertainty in 
the selection of host halo masses. For the Menc i W DM model, the shaded region also accounts for the unce rtainties in the suppres sion of 
progenitors below the free-streaming scale (see I3.2( l The points with vertical error bars are measurements (Nierenberg et al. 2012). The 
top panels show the comparison for lower mass hosts at lower (left) and higher redshift (right), while the lower panels are the equivalent 
for higher mass hosts. 



In the previous section we found significant differences 
in the model predictions for the satellite luminosity func- 
tion. As discussed in Section [31 numerous physical pro- 
cesses contribute to the final predicted luminosity func- 
tion. In this section, we compare key aspects of the mod- 
els in order to elucidate which model assumptions drive 
the predicted differences. We first compare the subhalo 
mass function in the models, before the effects of tidal 
stripping take effect, for all subhalos regardless of final lu- 
minosity. We then show the stellar mass to halo mass re- 
lationship for hosts and satellites in the luminosity range 
considered in this work, before the subhalos have un- 
dergone tidal stripping. We conclude by comparing the 
predicted satellite colors to show the differences in star 
formation histories. 

5.1. Unstripped subhalo mass function 

The four simulations all assumed slightly different cos- 
mologies and had different mass resolutions, so it is im- 
portant to consider how much the differences in the pre- 
dicted results are due to these effects. In Figure[3jwe plot 
the mass function of subhalos in a mass range represen- 
tative of the hosts and subhaloes from Figure [2j Note 



that this mass function is purely based on dark matter 
mass, and does not directly correspond to the satellite 
galaxies in Figure [21 as there is no luminosity selection. 

The Guo, Menci, and Lu CDM subhalo mass functions 
are extremely similar, while the WDM model shows sig- 
nificant truncation for subhalos with \og[M su b/Mhost\ < 

— 1 for the less massive hosts and log[M s „h/M^ ost ] < 

— 1.5 for the more massive hosts, and predicts similar 
behavior to other WDM models with sub-Mpc cutoff 
scales (e.g. ICoh'n et all l2000t ISmith fc Markovid l201ll : 
iKamada et alJl2013t ). 

5.2. The halo mass to stellar mass relationship 

One of the most important physical processes relied on 
in this comparison between simulation and observation, 
is the relationship between host halo mass and stellar 
mass, because the number of subhalos around a given 
host galaxy is strongly dependent on the halo mass of 
the host galaxy (see e.g. lBusha et al.ll2011[ ). particularly 
in the case of CDM. Furthermore, as we always consider 
the quantity log 10 [L s /L^] > —3, it is important to dis- 
tinguish whether differences in the models are caused by 
differences in L s or in L^. 
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Fig. 3. — Comparison of the subhalo mass function at the time of subhalo accretion for the four models, for a range of masses representative 
of the host halos and subhalos in Figure [2] Note that these mass functions do not include the effects of tidal stripping by the main halo, 
nor are they selected to host luminous satellite galaxies, or host galaxies with the same stellar masses as the hosts in Figure [2] 



In Figure [4] we plot the halo to stellar mass relation- 
ship for the four models, in addition to the observed 
relationship from iDutton et al.l (|2010() , which was used 
to estimate R200 for the host galaxies in the observa- 
tions. For bright galaxies with stellar masses greater than 
1O 1O ' 5 M0, the models are all very similar to each other 
and to the observed relationship, within the large scat- 
ter. Thus we conclude that the differences in amplitude 
in the predictions for the satellite luminosity function are 
not driven primarily by differences in the halo to stellar 
mass relationship for host galaxies. 

Towards the faint end, all models are very similar 
within the scatter, although the Menci models show a 
marginally higher dark matter mass at fixed stellar mass. 
As we show in the following section, the Menci satellite 
colors tend to be bluer than the Guo and Lu colors. From 
(|Bell et al.ll2003ft . Table 7, the bluer colors of the Menci 
model, correspond to an r band mass-to-light ratio which 
is lower by roughly 0.1-0.2 dex on average than for the 
Guo and Lu models. Thus, as the final luminosity func- 
tions are similar, it is expected that the Menci models 
would produce galaxies with on average slightly lower 
stellar masses for fixed halo mass. 



5.3. Satellite galaxy colors 



The intrinsic colors of satellite galaxies are dependent 
on the star formation history. As discussed above, this 
is determined by the metallicity, feedback and UV heat- 
ing, in addition to environmental effects as the satellites 
enter the influence of the host galaxy halo. All three 
models u sed the same stellar popu lation synthesis mod- 
els from ([Bruzual fc Ch ariot 2003|), thus satellite galaxy 
colors provide an important means of distinguishing be- 
tween different physical models for the suppression of 
star formation in low mass halos. Two models which 
produce similar luminosity functions may rely on very 
different star formation prescriptions, which will result 
in different color distributions. In Figure [SJ we plot the 
predicted distribution of rest-frame u-i colors for satel- 
lites with log 10 [L s /Lh] > —3. 

In order to summarize these predictions, we fit each 
of the distribu tions to a skewed normal distribution 
(|Azzalinilll985l ) where the probability of having color c is 
given by: 

P(c) «e^(l + er/(f^)) (1) 

Here c Q and a are the usual mean and standard de- 
viations of a normal distribution while the parameter a 
describes the skewness. The best fit values for these pa- 
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Fig. 4. — Comparison of the virial to stellar mass relationship for both host and satellite galaxies for the four models. The subhalo virial 
masses are taken at the time of accretion, and thus do not include the effects of tidal stripping. 



rameters are listed in Table Q] 6 . The parametrization of 
the distributions as being skewed normal is not physi- 
cally motivated, but rather intended to facilitate future 
comparison between observation and these predictions. 

The models all predict significantly different color dis- 
tributions for the satellite galaxies, with the exception of 
the Menci CDM and WDM models, for which the same 
star formation parameters were used. We highlight the 
fact that although the luminosity function predicted in 
the Menci WDM and Lu CDM models are similar within 
some redshift and host stellar mass ranges, the color dis- 
tributions are very different. This is due to the fact that 
in the Lu CDM model, the faint end of the luminosity 
function is suppressed mostly by the effects of feedback 
and ram pressure stripping and heating of the gas by 
the host halo, while in the Menci WDM model, the lu- 
minosity function is suppressed by the lack of low mass 
subhalos, as expected. 

In order to facilitate comparison with observation, it is 
important to also consider the effects of dust extinction. 
In Figure [6j we show the prediction for the rest-frame 
colors of satellite galaxies with dust. The Menci models 
behave very differently with the addition of dust, relative 
to the Guo and Lu models, with the Menci color distri- 

6 We do not consider the secondary peaks in the Menci models 
in this fit. 



butions becoming much wider and significantly redder 
while the Guo and Lu models develop a longer redward 
tail, without significant other alteration. The secondary 
peaks in the Menci model become more prominent with 
the addition of dust, so we fit the Menci models with the 
sum of two skewed normal distributions with parameters 
reported in Table [U where the value A2/A1 describes the 
relative amplitudes between the two skewed normals. 

6. COMPARISON WITH LITERATURE 

The models presented in this work predict a broad 
range of observables, many of which have already been 
compared to observations, either at lower redshifts, or 
of brighter objects than those considered in this pa- 
per. Below we provide a brief comparison with some of 
these observations, a more detailed discussion of many of 
these comparisons c an be found in the original papers de- 
scribing the m odels (jMenci et al.ll2012tlGuo et al.ll2011bt 
ILu et al.ll20ll . 

6.1. Satellite luminosity function 

Studies of the satellite luminosity function typically 
focus on satellite galaxies at low redshifts with z < 0.1, 
bo th in simulation and obs ervations. As we have shown 
in iNierenberg et~all (120121 1. the COSMOS field can be 
used to study satellites at higher redshifts, and yields 
luminosity functions at low redshifts which are con- 
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TABLE 1 

Best fit parameters of the color distributions 



Host Properties 


Dust 


Guo et al. 2012 


Lu et al. 2012 


Menci, CDM 


Menci, WDM 






c , cr, a 


c , cr, a 


c , cr, a, 

A2/A1, c , 2, o" 2 , a 2 


c , cr, a, 

A2/A1, c , 2, o"2, a 2 


10.5 < Iog[M*/M ]< 11.0, 0.1<z<0.4 


Yes 


1.7, 0.5, 8.8 


1.9, 0.4, 2.5 


1.6, 0.2, -5.0 
1.4, 1.0, 1.5, 12.5 


1.4, 1.0, 9.3 

0.2, 3.3, 0.25, -1.9 




No 


2.0, 0.3, 1.5 


1.7, 0.5, 7.0 


1.3, 0.3, 1.4 


1.9, 0.5, -4.0 


10.5 < log[M*/M©]< 11.0, 0.4<z<0.8 




1.6, 0.5, 5.6 


1.8, 0.3, 2.4 


1.4, 0.3. -9.0 
3.1, 2.9, 0.9, 2.5 


1.3, 0.3, 3.0 
1.5, 2.1, 0.5, 1.6 






1.9, 0.3, 1.2 


1.6, 0.4, 4.5 


1.7, 0.5,-5.4 


1.7, 0.4, -4.4 


11.0 < Iog[M*/M ]< 11.5, 0.1<z<0.4 




1.8, 0.6, 5.1 


1.2, 0.5, 2.7 


1.5, 0.24, -5.9 
4, 2.6, 0.5, 


1.3, 1.1, 5.5 
0.3, 3.3, 0.3,-4.5 






2.5, 0.3, -1.2 


1.8, 0.5, 4.3 


2.0, 0.5,-4.4 


1.9, 0.5, -4.2 


11.0 < \og[M*/M Q }< 11.5, 0.4<z<0.8 




1.7, 0.5, 4.6 


1.9, 0.5, 2.9 


1.3, 0.25, 

2.3, 2.9, 0.7, -2.6 


1.6, 0.2, 

1.3, 2.0, 0.5, 1.7 






2.4, 0.3, -1.3 


1.7, 0.5, 3.7 


1.7, 0.5,-6.2 


1.7, 0.4, -5.3 



10 



Nierenberg et al. 



2.5 

2.0 

1.5 

1.0- 

0.5- 



1 1 1 1 1 1 r 

10.5<log 10 [M;/Mo]<11.0 



0.1 <z 



a o.o 

o 2.5 

u 
as 



1.5- 
1.0- 
0.5- 
0.0 




0.4<z<0.8 




ll.O<log lo [M ft 7M0]<11.5 


1 1 1 1 1 1 1 

Lu et al. 2012 
■M Guo et al. 2011 
" Mencl, CDM 

<MM Menci, WDM 


. ' - - -^ -i 


s 1 1 1 — .1 - f — 





0.5 1.0 1.5 2.0 2.5 3.0 3.5 0.5 1.0 1.5 2.0 2.5 3.0 3.5 

u-i (Rest, with dust) 

Fig. 6. — The normalized rest-frame distribution of colors, with dust extinction for the satellite galaxies plotted in Figure 2, along with 
the best fitting skewed normal distributions. The Menci models are fit with a sum of two skew normal distributions, to account for the 
prominent secondary peak that appears. 
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sistcn t with satellite luminosity functions from SPSS 
(e.g. IGuo et al.ll2011al ILlu et alJl20ll ILares et alJl20ll 
iStrigari fc Wechslerll2012D. and the Milky Way satellite 
luminosity function ()Tollerud et al.l 120081) . All four of 
the models provide good fits to the satellite luminosity 
function of Milky Way mass hosts at low redshifts, by de- 
sign. The extra information of redshift and host galaxy 
stellar mass provides additional constraints to these mod- 
els. For instanc e, when comparing to bright f ield galax- 
ies obse rved by iPerez-Gonzalez et al.l (|2008l) , IGuo et ah! 
(|2011b|) found their model agreed well at low redshift, 
while becoming discrepant with observation by a redshift 
ofl. 

6.2. Satellite Colors 

In Section 15.31 we showed that the three CDM mod- 
els predicted significantly different distributions for the 
colors of faint satellite galaxies. Although we cannot di- 
rectly test this prediction in this work, some comparisons 
can be made with measu rements of field ga laxies at low 
redshifts. In particular, IGuo et aTl (|2011b| ) found that 
their model predicted colors matched observed SDSS col- 
ors well except for low masses (log 10 [M*/Mq] < 9.5), 
which were redder in th e simulation than in observa- 
tions. M enci et al.l <|201 2h found that down to an absolute 
magnitude of M r = — 18 (roughly a stellar mass of ^9 
log[A/*/M©]), their color distribution agreed well with 
SDSS measurements from lBaldrv et al.l (|2004f >. They did 
not compare for fainter satellites. 

The Lu and Guo semi-analytic models were not opti- 
mized to reproduce the observed galaxy color-magnitude 
relation at low redshifts, thus future implementations of 
these semi-analytic models which include this informa- 
tion may yield significantly different predictions for the 
color distri bution of satel li te gal axies. 

Recently, iKnobel et al.l (|2012l ) measured the red frac- 
tion of massive (log 10 [M*/Mq] > 10) satellites between 
redshift 0.1 and 0.8, and found no evidence for signifi- 
cant evolution over this time, indicating that these low 
redshift results may apply to higher redshifts. Future 
measurements of faint satellite colors since intermedi- 
ate redshifts will provide interesting new constraints on 
these models, by helping to distinguish between suppres- 
sion of the satellite luminosity function by environmental 
quenching and supernovae feedback, or the subhalo mass 
function at the low mass end. 

7. DISCUSSION AND CONCLUSION 

We have compared predictions from three current semi- 
analytic models applied to CDM simulations and one 
semi-analytic model applied to a WDM simulation to 
observations of satellite galaxies in two redshift intervals 
and with two host mass bins. By comparing multiple 
models, we demonstrated that the effects of varying star- 
formation prescriptions within semi-analytic models, and 
varying the underlying dark matter mass function have 
fundamentally different effects on the predicted host stel- 
lar mass dependence, and redshift evolution of the satel- 
lite luminosity function. In particular, we show that the 
WDM model predicts a satellite luminosity function with 
much weaker host stellar mass dependence and redshift 
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evolution than any of the CDM models. 

We find that the WDM provides the best match to 
observation in all redshift and host mass intervals, most 
closely matching the host mass dependence and lack of 
redshift evolution in the data. This comparison high- 
lights the importance of comparing models of satellite 
galaxy evolution to observations from a range of red- 
shifts and stellar mass regimes, as a model that provides 
a close match to observation in one regime can perform 
more poorly than other models in a different regime. 

This exploration of different CDM models suggests 
that current models generically have difficulty reproduc- 
ing the mass and redshift dependence of the satellite 
luminosity function. However, future improvements to 
the semi-analytic models for star formation prescriptions 
may allow for improvements between the CDM predic- 
tions of the satellite l uminosity function. In fact, as 
shown in the paper by iNierenberg et al.l (l2012f) sub halo 
abundance matching techniques (|Busha et al.ll2011| ) can 
reconcile our observations with the CDM subhalo mass 
function. Although abundance matching is descriptive 
and not directly linked to known physical processes, it is 
possible that future semi-analytic models will be able to 
produce the stellar mass to halo mass relation as illus- 
trated by this approach. 

Additional observational data are needed to further 
constrain the models and therefore help distinguish be- 
tween improved baryonic physics in CDM models and 
a WDM power power spectrum. For example, we have 
shown in this work that the Lu CDM model which most 
closely matches the data predicts a significantly differ- 
ent color distribution for the satellite galaxies than the 
Menci WDM model, indicating that future observations 
of the colors of faint satellites will provide an impor- 
tant test of whether the luminosity function of satellites 
is suppressed primarily by baryonic processes or by a 
WDM mass function. We plan to carry out such mea- 
surement in the near future by exploiting the rich multi- 
color datasets publicly available in the HST legacy fields. 
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